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Abstract. We study paramctrization of polarized parton distribution functions in the 
a s leading order (LO) and in the next-to-leading order (NLO). From x 2 fitting to the 
experimental data on A\, optimum polarized distribution functions are determined. 
The quark spin content AE is very sensitive to the small- a; behavior of antiquark 
distributions which suggests that small- a; data are needed for precise determination of 
AS. We propose three sets of distributions and also provide FORTRAN library for 
our distributions. 



INTRODUCTION 

Experimental data on polarized structure functions g\ have been accumulated 
for the last several years. The data with the proton, deuteron, and 3 He targets are 
now available and new data are expected to be given by RHIC at BNL, COMPASS 
at CERN, and other facilities in the near future. In the light of such progresses 
and future projects, we should summarize present knowledge of polarized parton 
distribution functions by using all available data at this stage. For this purpose, 
we formed the group called Asymmetry Analysis Collaboration (AAC) and tried 
to determine the polarized parton distribution functions. 

In this paper, we explain our analysis to determine optimum polarized distribu- 
tions by x 2 fitting to the experimental data of spin asymmetry A\. In addition to 
the discussions based on the work in Ref. [1], we introduce a FORTRAN library 
for our distributions as a recent progress. 



1 ) 98td25@edu.cc.saga-u.ac.jp 

2 ) hyuki@bnl.gov 

3 ) miyama@comp.metro-u.ac.jp 



PARAMETRIZATION 



We determine initial polarized parton distributions at Q 2 = 1.0 GeV 2 (= Qq). 
Considering the counting rule, we adopt the following functional form: 

Afi(x, Ql) = A, x a > (1 + 7l x A fi(x, Ql). (1) 

Here, Afi and represent the polarized and unpolarized parton distributions, 
respectively and A iy a iy 7,, and A, are free parameters. These parameters are 
constrained by the positivity condition, | Afi(x) | < fi(x). We simply apply it 
not only for the leading-order (LO) case but also for the next-to-leading-order 
(NLO) case although it is valid only in LO. Furthermore, we assume SU(3) flavor- 
symmetric sea because we don't have enough data to extract flavor dependence of 
polarized sea. Under this assumption, the first moments of Au v and Ad v can be 
fixed by the axial charges for octet baryon, F and D. Therefore, we should determine 
Au v , Ad v , Aq, and Ag distributions and total number of the parameters becomes 
14. 

As the experimental data to which the parameters are fitted, we chose spin 
asymmetry A\ instead of g± since A\ is closer to direct observable in experiment 
rather than g±. Ai is given by 

A, (x Q 2 ) ~ 91{X ' Q2) - a, (x C 2- M^ + R(*,Q 2 )] (2) 

where, R is the ratio of longitudinal to transverse cross-sections. In our analysis, 
Ai is calculated by using F 2 obtained by GRV98 distributions, g\ obtained by 
our parametrized distributions, and R taken from the SLAC-1990 analysis. Then, 
X 2 = E[Ai^(x,Q 2 ) - A'f c (x,Q 2 )} 2 /[<j data (x,Q 2 )} 2 is calculated and is minimized 
by the subroutine MINUIT. Here, Af ata , Al alc , and <j data indicate the experimental 
Ai, calculated A 1 , and experimental error, respectively. 



RESULTS AND DISCUSSIONS 

From the \ 2 analysis, we obtain the results with x 2 = 322.6 in the LO case, 
and x 2 = 300.4 in the NLO case for 375 data points. The NLO x 2 is significantly 
smaller than the LO one. It suggests that the NLO analysis is necessary for precise 
analysis. The obtained A 1 for the proton and for the neutron at Q 2 = 5.0 GeV 2 
are shown in Fig.l with the experimental data. Although the comparison is not 
straightforward because the data are taken at various Q 2 , the obtained parame- 
ters reproduce the data well both in the LO and the NLO case. Figure 2 show 
the obtained polarized parton distribution functions at Q 2 = 1.0 GeV 2 . In our 
analysis, the first moment of Au v (Ad v ) is fixed by positive (negative) value and 
the obtained distribution becomes positive (negative). Furthermore, antiquark and 
gluon distributions become negative and positive, respectivly. 
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FIGURE 1. Spin asymmetries A\ for the proton (left figure) and neutron (right figure). 

For these distributions, quark spin content becomes AS = 0.201 for LO and 
AS = 0.051 for NLO at Q 2 =1.0 GeV 2 . The NLO AS is significantly smaller than 
the LO one. It is also smaller than the generally quoted values in other analyses that 
range in 0.1 ~ 0.3. The difference of AS comes mainly from the difference of Aq in 
the small- a; region. Because we don't have a small- a; data at this stage, we should 
fix the small- a; behavior of Aq from theoretical suggestions. We perform x 2 fitting 
by fixing the parameter aq, which controls the small-x behavior of Aq, to aq = 1.0 
and 1.6 according to the predictions of the Regge theory and the perturbative QCD. 
These values are rather larger than the value aq = 0.32 ± 0.22 which is obtained 
by the previous NLO analysis. As a result, we get AS = 0.241 (x 2 = 305.8) for 
a q = 1.0 and AS = 0.276 (x 2 = 323.5) for a- q = 1.6. The obtained AS are larger 
than the previous NLO result and are within the usually quoted range. From these 
results, we find that AS is very sensitive to the small-a; behavior of the antiquark 
distributions. In order to determine the quark spin content precisely, the data in 
the small- a; region are needed. Because x 2 f° r the aq = 1.0 case is close to the one 
for the previous NLO, this set can be also considered as a good parametrization. 
Therefore, we propose three sets of distributions, LO, NLO with free aq (NLO-1), 
and NLO with aq = 1.0 (NLO-2), as AAC polarized distribution functions. 

Finally, we briefly comment on the error estimate of the obtained distributions 
[2]. The uncertainties of our distributions can be estimated by using the error 




FIGURE 2. Obtained polarized distributions at Q 2 = 1.0 GeV 2 . 
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matrix which is obtained by MINUIT. The 1-a boundary of the distribution F is 
given by 




(3) 



where, Vij is the error matrix element for parameters and aj. This analysis is in 
progress and will be reported elsewhere. 



For practical application of our distributions, we provide a library program which 
is given as FORTRAN subroutine AACPDF(ISET, Q2, X, POLPDF, STRUCT). 
The users should call this subroutine in their programs with the input parameters 
ISET, Q2, and X. This library contains three sets of distributions, LO, NLO-1, and 
NLO-2, and the parameter ISET designates which set is used. When ISET=1, 2, 
or 3, the returned values are those in the LO, NLO-1, or NLO-2, respectively. The 
parameters Q2 and X specify Q 2 and Bjorken-rr at which the distributions are calcu- 
lated. The allowed ranges are 1.0 GeV 2 < Q 2 < 10 6 GeV 2 and 1(T 9 < x < 1.0. The 
obtained values are returned by the arrays POLPDF(-3:3) and STRUCT(3). The 
available polarized parton distribution functions and the structure functions are 
listed in Table 1. It should be noted that the returned values are the distributions 
and the structure functions multiplied by x. 

The distribution values are obtained by interpolating the grid data which are 
provided as DATA statements in the subroutine. The number of the grid points 
is 23 for the variable Q 2 and 68 for the variable x. Because the Q 2 dependence of 
the distributions is almost linear function of t = InQ 2 , we simply use the linear 
interpolation of t for the variable Q 2 . If t' is in the range U < t' < t i+1 , the 
distributions at t — t' are approximated by 



LIBRARY 



f(t') 




f(U) + 



f-U 



f(t i+ i). 



(4) 



TABLE 1. Available distributions and structure functions. 



I 



POLPDF(I) 



STRUCT (I) 



-3 
-2 
-1 

1 
2 
3 



s quark [xAs(x, Q 2 ) = xAs(x, Q 2 ))] 

d quark [xAd(x, Q 2 ) — xAd sea (x, Q 2 )] 

u quark [xAu(x,Q 2 ) = xAu sea (x,Q 2 )] 

gluon [xAg(x,Q 2 )} 

u- valence quark [xAu„(x, Q 2 )] 

d- valence quark [xAd v (x, Q 2 )] 

s quark \xAs(x,Q 2 ) = xAs(x, Q 2 )] 



xgi(x,Q 2 ) for proton 
xgi(x,Q 2 ) for neutron 
xgi(x, Q 2 ) for deutcron 



4 



On the other hand, the x dependence of the distributions is rather complicated. 
Therefore, we use the cubic spline interpolation for the variable x. If x' is in the 
range xt < x' < x i+ i, the distributions at x — x' are approximated by 

f(x') « f(xi) + Bi(x' - Xi ) + Ci(x' - Xl f + Di(j> - Xl f. (5) 

Here, Bi, Cj, Di are the spline coefficients which can be calculated by the grid data. 
By using these interpolation methods and the grid data, the AACPDF returns the 
values of the polarized parton distribution functions and the structure functions g± 
at the specified Q 2 and x point. We verified that the distributions obtained by the 
library reproduce well the results in Ref. [1]. 

AAC library can be downloaded from our web page [3]. The subroutine AACPDF 
is in the file " aac.f " . This file also includes the subroutine SPLINE and the function 
ISERCH which are used in the AACPDF. 

SUMMARY 

The polarized parton distribution functions have been determined by \ 2 fitting 
to Ai experimental data. As a result, we found that the NLO x 2 is significantly 
smaller than the LO one. It implies that the NLO analysis is important. In the 
NLO analysis, the obtained AE is rather smaller than the usually quoted values, 
and the differences come from the small-x behavior of Aq. In order to fix AS 
precisely, the small- x measurements are required. We propose three sets of AAC 
distributions, LO, NLO-1 (a q : free), and NLO-2 (a q = 1.0 fixed). The FORTRAN 
library is provided for calculating the AAC distributions numerically. This library 
is useful for practical application of the AAC distributions and available at our web 
site [3]. 
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